Spectrum Sensing in Cognitive Radios Based on 
Multiple Cyclic Frequencies 



Invited Paper 



o 
o 

(N 



C/5 

O 



> 
&\ 
O 
On 
O 

o 
o 



X 



Jarmo Lunden*, Visa Koivunen*^, Anu Huttunen* and H. Vincent Poor^ 
* SMARAD CoE, Signal Processing Laboratory 
Helsinki Univ. of Technology, Finland 
Email: jrlunden@wooster.hut.fi 
^School of Engineering and Applied Science, 
Princeton University 
Email: poor@princeton.edu 



Abstract — Cognitive radios sense the radio spectrum in order 
to find unused frequency bands and use them in an agile manner. 
Transmission by the primary user must be detected reliably even 
in the low signal-to-noise ratio (SNR) regime and in the face 
of shadowing and fading. Communication signals are typically 
cyclostationary, and have many periodic statistical properties 
related to the symbol rate, the coding and modulation schemes 
as well as the guard periods, for example. These properties 
can be exploited in designing a detector, and for distinguishing 
between the primary and secondary users' signals. In this paper, 
a generalized likelihood ratio test (GLRT) for detecting the 
presence of cyclostationarity using multiple cyclic frequencies is 
proposed. Distributed decision making is employed by combining 
the quantized local test statistics from many secondary users. 
User cooperation allows for mitigating the effects of shadowing 
and provides a larger footprint for the cognitive radio system. 
Simulation examples demonstrate the resulting performance 
gains in the low SNR regime and the benefits of cooperative 
detection. 

I. Introduction 

Spectrum sensing is needed in cognitive radios in order 
to find opportunities for agile use of spectrum. Moreover, 
it is crucial for managing the level of interference caused 
to primary users (PUs) of the spectrum. Sensing provides 
awareness of the radio operating environment. A cognitive 
radio may then adapt its parameters such as carrier frequency, 
power and waveforms dynamically in order to provide the best 
available connection and to meet the user's needs within the 
constraints on interference. 

In wireless communication systems we typically have some 
knowledge on the waveforms and structural or statistical 
properties of the signals that the primary user of the spectrum 
is using. Such knowledge may be related to the modulation 
scheme, the symbol or chip rate of the signal, the channel 
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coding scheme, training or pilot signals, guard periods, and 
the power level or correlation properties of the signal, just 
to mention a few. These properties may be used to design a 
detector that works in a very low SNR regime and has low 
complexity and consequently low power consumption. These 
are very desirable properties especially for cognitive radios 
in mobile applications. In the absence of any knowledge of 
the signal, one may have to resort to classical techniques 
such as energy detection [ 1 ] . An energy detector may need to 
collect data over a long period of time to detect the primary 
users reliably. Moreover, controlling the false alarm rates in 
mobile applications is difficult because the statistics of the 
signals, noise and interference may be time-varying. Another 
significant drawback is that energy detection has no capability 
to distinguish among different types of transmissions or to 
dichotomize between primary and secondary users of the 
spectrum. 

Cyclostationary processes are random processes for which 
the statistical properties such as the mean and autocorrelation 
change periodically as functions of time [2]. Many of the 
signals used in wireless communication and radar systems pos- 
sess this property. Cyclostationarity may be caused by modu- 
lation or coding |2|, or it may be also intentionally produced 
in order to aid channel estimation or synchronization Q. 
Cyclostationarity property has been widely used in intercept 
receivers |2|, [4|, |5|, direction of arrival or time-delay estima- 
tion, blind equalization and channel estimation 1 6 1 as well as in 
precoder design in multicarrier communications J3). In order 
to exploit cyclic statistics, the signal must be oversampled with 
respect to the symbol rate, or multiple receivers must be used 
to observe the signal. The use of cyclostationary statistics is 
appealing in many ways: noise is rarely cyclostationary and 
second-order cyclostationary statistics retain also the phase 
information. Hence, procedures based on cyclostationarity tend 
to have particularly good performance at the low SNR regime. 
Moreover, cyclostationarity allows for distinguishing among 
different transmission types and users if their signals have 
distinct cyclic frequencies. A comprehensive list of references 



on cyclostationarity along with a survey of the literature is 
presented in Q. 

The presence of cyclostationary signals may be determined 
by using hypothesis testing. Many existing tests, such as 
|8l , are able to detect the presence of cyclostationarity at 
only one cyclic frequency at a time, and they partly ignore 
the rich information present in the signals. For example, a 
communication signal may have cyclic frequencies related to 
the carrier frequency, the symbol rate and its harmonics, the 
chip rate, guard period, the scrambling code period, and the 
channel coding scheme. In this paper we propose a method 
for detecting multiple cyclic frequencies simultaneously. It 
extends the method of |8| to take into account the rich 
information present at different cyclic frequencies. This pro- 
vides improved detector performance over techniques relying 
only on single cyclic frequency and facilitates dichotomizing 
among the primary and secondary user signals and different 
waveforms used. 

In cognitive radio systems, there are typically multiple geo- 
graphically distributed secondary users (SUs) that need to de- 
tect if the primary user is transmitting. The distributed sensors 
may work collaboratively to decide between two hypotheses: 
is the primary user active, or is the spectrum unused and 
available for the secondary users? Decentralized processing 
has a number of advantages for such situations. Obviously, 
it allows for a larger coverage area. Furthermore, there are 
gains similar to diversity gains in wireless communications 
so that the detection becomes less sensitive to demanding 
propagation conditions such as shadowing by large obstacles, 
large numbers of scatterers, differences in attenuation, or 
fast fading caused by mobility. Moreover, distributed sensory 
systems may require less communication bandwidth, consume 
less power, be more reliable and cost less as well. In this 
paper, we propose a simple decentralized decision making 
approach based on sharing and combining quantized local 
decision statistics. This approach may be used in both decision 
making with or without a fusion center. 

This paper is organized as follows. In Section II, there is a 
short review of cyclostationary statistics. A novel detector for 
multiple cyclic frequencies is derived in Section III. Section 
IV addresses the problem of collaborative detection of primary 
user. Simulation results demonstrating the detector's reliability 
in the low SNR regime as well as the gains obtained via 
collaborative operation are presented in Section V. Finally, 
conclusions are drawn in Section VI. 

II. Cyclostationarity: a recap 

In this section, we provide a brief overview of cyclostation- 
arity in order to make the derivation of the detector in Section 
III clearer. A continuous-time random process x(t) is wide 
sense second-order cyclostationary if there exists a Tq > 
such that 0: 

=l*x(t + T ) Vf (1) 

and 

R x {tiM = R x (t 1 +T ,t 2 +T ) Vti,t 2 . (2) 



To is called the period of the cyclostationary process. 

Due to the periodicity of the autocorrelation R x (ti,t2), it 
has a Fourier-series representation. By denoting t\ — t + r/2 
and ti =t — r/2, we obtain the following expression for the 
Fourier-series 0: 

RAt+ T ^t- T -)=Y J K{r)e^ a \ (3) 

a 

where the Fourier coefficients are 

1 r°° t t 

and a is called the cyclic frequency. The function R x (t) is 
called the cyclic autocorrelation function. If the process has 
zero mean, then this is also the cyclic autocovariance function. 

When the autocorrelation function has exactly one period 
T we have the following set of cyclic frequencies 

A= {a = k/T,k > 1}, 

where R x (t) is the cyclic autocorrelation function and A 
are the set of cyclic frequencies. The cyclic frequencies are 
harmonics of the fundamental frequency. If the autocorrelation 
function has several periods To, Xi, . . ., we may express R"(t) 
at the limit 

R a x {r) = lim i f T/2 x(t + T -)x*{t - T -)e^dt. (5) 

T^oo 1 J_ T /2 I I 

The process x(t) is almost cyclostationary in the wide sense 
and the set of cyclic frequencies A is comprised of a countable 
number of frequencies that do not need to be harmonics of the 
fundamental frequency. In general, the process is said to be 
cyclostationary if there exists an a ^ such that R"(r) ^ 
for some value of r. Typically cyclic frequencies are assumed 
to be known or may be estimated reliably. 

III. Detection using multiple cyclic frequencies 

Statistical tests for the presence of a single cyclic frequency 
have been proposed, for example, in (8). The tests in [8 | have 
asymptotically constant false alarm rate (CFAR) for testing 
presence of cyclostationarity at a given cyclic frequency. 
However, the tests do not retain the CFAR property over a 
set of tested frequencies. 

Typical communication signals exhibit cyclostationarity at 
multiple cyclic frequencies instead of just a single cyclic 
frequency. That is, for example a signal that is cyclostationary 
at the symbol frequency is typically cyclostationary at all 
integer multiples of the symbol frequency as well. There also 
may be cyclic frequencies related to the coding and guard 
periods, or adaptive modulation and coding may be used. In 
such cases the cyclic frequencies present may vary depending 
on channel quality and the waveform used. If one is testing 
for the presence of many different signals at a given frequency 
band, or in case the cyclic frequencies are not known, it would 
be desirable to retain the CFAR property over the whole set of 
tested cyclic frequencies. This would be especially desirable 
in a cognitive radio application where the interest is in finding 



unoccupied frequency bands. Otherwise the frequency band 
may unnecessarily be classified as occupied for most of the 
time. 

In the following we extend the test based on second-order 
cyclic statistics of [8| to multiple cyclic frequencies. To do so 
we first define all the terms used in the test statistics. 

Let (*) denote an optional complex conjugation. The no- 
tation allows convenient handling of both cyclic autocorre- 
lation and conjugate cyclic autocorrelation with only one 
equation. An estimate of the (conjugate) cyclic autocorrelation 
R xx (*) (a, t) may be obtained using M observations as 

M 



RxxC) 0, T ) = ^2 x(t)xM (t + r)e 
= ^w(a,T) +e(a,r), 



(7) 



where the latter term is the estimation error. This estimator is 
consistent, (see ID) so that the error goes to zero as M — ► oo. 

Now we need to construct a test for a number of lags 
T\, . . . , rjv as well as a set of cyclic frequencies of interest. 
Let A denote the set of cyclic frequencies of interest, and 

r xx (,)(a) = Re{E xa . ( *)(a,Ti)}, . . . , Re{R xx( ,) (a, t n )}, 

lm{R xxM (a,n)}, . . .,lm{R xxM (a,T N )} 

(8) 

denote a 1 x 2N vector containing the real and imaginary parts 
of the estimated cyclic autocorrelations at the cyclic frequency 
of interest stacked in a single vector. 

The 2N x 2N covariance matrix of r xx (») can be computed 
as © 
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(9) 



where the (m, n)th entries of the two covariance matrices Q 
and Q* are given by 



Q{m,n) = Sf Tm f Tn (2a,a) 



and 



(10) 



Q*(m,n) = Sl mfTn (0,-a). 

Here, Sf Tm f Tn (a,ui) and Sj r t (a,oj) denote the unconju- 
gated and conjugated cyclic spectra of f(t, t) = x{t)x^ (t + 
r), respectively. These spectra can be estimated using fre- 
quency smoothed cyclic periodograms as 

1 (L-l)/2 

— y w( s ) 
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where F T {io) = Ylf=i x{t)x& (t + T)e~ lwt and W is a 
normalized spectral window of odd length L. 

Now the hypothesis testing problem for testing if a is a 
cyclic frequency can be formulated as [8| 

H Q : V{r„}^ =1 => r xxM (a) = e xxM (a) 



Hi : for some {r n }^ =1 
=> f xx <,) (a) = r x 



(13) 



.(,)(«) +e xxM (a). 



Here e xx (») is the estimation error which is asymptotically nor- 
mal distributed, i.e., limM->oo V^e xx (») — N(0, X xx (*)) 0. 
Hence, using the asymptotic normality of r xx (*> the general- 
ized likelihood ratio (GLR) is given by 

-1 ~T 

A 



cxp( 



xx ( *> xxy* 



exp(-|M(r xx( „) - r xx (,))£ xxW (r x:r c,) - r^.)) 71 ) 
= exp^-Mr^wS-^,)^*))- 

(14) 

Finally, by taking the logarithm and multiplying the result 
by 2, we arrive at the test statistic in (SI 

T xxM (a) = -2 In A - Mf xxM t-^ xM . (15) 

Under the null hypothesis T xx (») (a) is asymptotically x\n 
distributed. 

Now in order to extend the test for the presence of second- 
order cyclostationarity at any of the cyclic frequencies of 
interest a £ A simultaneously, we formulate the hypothesis 
testing as follows 

Ho : Va £ A and V{T n } n=1 r xx (*) (a) — e xx («o (a) 
Hi : for some a £ A and for some {7"n}n=i 
=> r xx( „) (a) = r Il( ,) (a) + e xxM (a). 

(16) 

For this detection problem, we propose the following two 
test statistics: 

V m = maxT a . lW (a) = maxMr xxM (a)t,-^, ) (a)fl x(t , ) (a) 



(17) 



Z? s = ^ 7; x (o (a) 



(18) 



M 



M 



The first test statistic calculates the maximum of the cy- 
clostationary GLRT statistic (fTBT l over the cyclic frequencies 
of interest A while the second calculates the sum. Assuming 
independence of cyclic autocorrelation estimates for different 
cyclic frequencies the test statistic T) s is the GLRT statistic. 
Depending on the signal and the set of tested cyclic frequen- 
cies the test statistics may have different performances. This 
requires further research. 

The asymptotic distribution of T> s is under the null hypoth- 
esis X2NN a w h ere N a is the number of cyclic frequencies in 
set A. This is due to the fact that the sum of independent chi- 
square random variables is also a chi-square random variable 
whose degrees of freedom is the sum of the degrees of freedom 
of the independent random variables. 



In the following we derive the asymptotic distribution of the 
test statistic T> m under the null hypothesis. As stated above, 
under the null hypothesis T xx m (a) is asymptotically x\n 
distributed. The cumulative distribution function of the chi- 
square distribution with 2N degrees of freedom is given by 



F(x,2N) 



l(N,x/2) 
T(N) 



(19) 



where ^f(k,x) is the lower incomplete gamma function and 
r(fc) is the ordinary gamma function. For a positive integer k 
the following identities hold: 



T(k) = (k- 1)! 



k—1 



7 



(fc, x )=r(fc)-(fc-l)!e-^- T . 



(20) 



(21) 



Hence, the cumulative distribution function of the chi-square 
distribution with 2N degrees of freedom is given by 



JV-l 



F(x,2N) = l-e- x/2 



{x/2) n 



(22) 



n=0 



The cumulative distribution function of the maximum of d 
independent and identically distributed random variables is the 
cumulative distribution function of the individual random vari- 
ables raised to the power d. Thus, the cumulative distribution 
function of the test statistic T> m is given by 

\ d 



N-l 



F Vm (x,2N,d) = l-e- x / 2 J2 



Or/2)* 



(23) 



n=0 



The corresponding probability density function is obtained 
by differentiating the cumulative distribution function, i.e., 



/^(*,2iM) = |(l-e-^£^ 



n=0 



n\ ^ (n-l)! 

\n=0 n=l v ; 



(24) 

Consequently, the null hypothesis is rejected if 
Fp m (T> m , 2N, N a ) > 1 — p where p is the false alarm 
rate and N a is the number of tested cyclic frequencies. 

IV. Cooperative detection 

User cooperation may be used to improve the performance 
and coverage in a cognitive radio network. The users may 
collaborate in finding unused spectrum and new opportunities. 
Many of the collaborative detection techniques stem from 
distributed detection theory; see iflOl , ifTTl . In cognitive radio 
systems, there are typically multiple geographically distributed 
secondary users that need to detect whether the primary user 
is active. All the secondary users may sense the entire band 
of interest, or monitor just a partial band to reduce power 
consumption. In the latter case each SU senses a certain part 
of the spectrum, and then shares the acquired information with 
other users or a fusion center. 



The cooperation may then be coordinated by a fusion center 
(FC), or it may take place in an ad-hoc manner without a 
dedicated fusion center. Here we assume that a fusion center 
collects information from all K secondary users and makes a 
decision about whether the spectrum is available or not. We 
assume that each secondary user sends a quantized version 
of its local decision statistics (such as the likelihood ratio) to 
the FC. In the case of very coarse quantization, binary local 
decision may be sent. To derive a test for the FC, we assume 
that the sensors are independent conditioned on whether the 
hypothesis Hq or Hi is true. Then the optimal fusion rule 
is the likelihood ratio test over the received local likelihood 
ratios lc 

K 

T K =Y[k. (25) 

i=l 

In case the secondary users send binary decisions, the sum 
of ones may calculated and compared to a threshold. Here, 
we consider the simplest way of making the decision using 
generalized likelihood ratios. Instead of using the product of 
the generalized likelihood ratios, we can employ the sum of 
generalized log-likelihood ratios. We propose the following 
test statistic for the hypothesis testing problem ( fT3T > 



K 



(26) 



and the following two for the hypothesis testing problem (Q~6| 



K 



D m k = max T^), (a) 

i=l 



(27) 
(28) 



u£A i=l 



where T^\ t) (a) is the cyclostationarity based test statistic (fT3T > 
from i*^ secondary user. Due to the use of generalized 
likelihood ratios, no optimality properties can be claimed. The 
GLRT test does, however, perform highly reliably in many 
applications. 

Under the conditional independence assumption the asymp- 
totic distributions of the test statistic and T> S: k are under 
the null hypothesis X2NK anc ^ X2NN a K> respectively. This 
is again due to the fact that the sum of independent chi- 
square random variables is also a chi-square random variable 
whose degrees of freedom is the sum of the degrees of 
freedom of the independent random variables. The cumulative 
distribution function of T> m ,K is under the null hypothesis 
F-p m (T) m k,2NK, N a ) where N a is again the number of 
tested cyclic frequencies. The testing is done similarly as in 
one secondary user case. 

Different techniques for reducing the amount of transmitted 
data, taking into account the relevance of the information 
provided by secondary users as well as how to deal with com- 
munication rate constraints will be addressed in a forthcoming 
paper. 



V. Simulation examples 

In this section the performance of the proposed detectors is 
considered. The test signal is an orthogonal frequency division 
multiplex (OFDM) signal. The baseband equivalent of a cyclic 
prefix OFDM signal may be expressed as 



= 0.05 



x(t) 



Nc-1 

E 

n=0 1= 



■IT.) 



(29) 



where N c is the number of subcarriers, T s is the symbol length, 
git) denotes the rectangular pulse of length T s , and the c n /s, 
denote the data symbols. The symbol length is the sum of 
the length of the useful symbol data T d and the length of the 
cyclic prefix T cp , i.e., T a = T d + T cp . 

The above OFDM signal exhibits cyclostationarity (i.e., 
complex conjugation is used in © and the following equa- 
tions) with cyclic frequencies of a = k/T s , k = 0, ±1, ±2, . . . 
and potentially other frequencies depending on the coding 
scheme. The cyclic autocorrelation surfaces for a = k/T s 
peak at r = ±T d (9). 

In the following the performance of cyclic detectors based 
on one and two cyclic frequencies is compared as a function of 
signal-to-noise ratio (SNR) in an additive white Gaussian noise 
(AWGN) channel. The SNR is defined as SNR = 10 log 10 
where o 2 x and cr^ are the variances of the signal and the noise, 
respectively. The cyclic frequencies employed by the detectors 
are 1 /T s and 2/T s . The detector based on one cyclic frequency 
uses the first frequency and the detectors based on two cyclic 
frequencies use both frequencies. Each detector uses two time 
lags ±T d . 

The cyclic spectrum estimates were calculated using a 
length-2049 Kaiser window with j3 parameter of 10. A Fast- 
Fourier transform (FFT) was employed for faster computation. 
The FFT size was 10000 giving a cyclic frequency resolution 
of 0.0001. 

The OFDM signal has 32 subcarriers and the length of the 
cyclic prefix is 1/4 of the useful symbol data. The subcarrier 
modulation employed is 16-QAM. The signal length is 100 
OFDM symbols. 

Fig- [U depicts the performance of the detectors as a function 
of the SNR for a constant false alarm rate of 0.05. Fig.|2]shows 
a zoom of the important area illustrating the differences in 
performance more clearly. All the curves are averages over 
10000 experiments. It can be seen that the detectors based on 
multiple cyclic frequencies outperform the detector based on 
single cyclic frequency in the low SNR regime. Furthermore, 
the multicycle detector calculating the sum over the cyclic 
statistics of different frequencies has the best performance. 

Fig. [3] plots the probability of detection vs. false alarm 
rate for SNR of -7 dB. The figure show that the detectors 
have desirable receiver operating characteristics. That is, the 
probability of detection increases as the false alarm rate 
parameter is increased. 

Next the performance gain from cooperative detection of 
several secondary users is analyzed. The signal is the same 
as above. The cooperative detection is based on the data of 5 
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Figure 1. Probability of detection vs. SNR. The multicycle detectors achieve 
better performance than the single cycle detector in the low SNR regime. The 
sum detector of the test statistic T> s has the best performance. 
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Figure 2. Probability of detection vs. SNR. Zoom of the important region. The 
multicycle detectors achieve better performance than the single cycle detector 
in the low SNR regime. The sum detector of the test statistic T> s has the best 
performance. 
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Figure 3. Probability of detection vs. false alarm rate. The detectors based 
on multiple cyclic frequencies achieve better performance than the detector 
based on a single cyclic frequency. 
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Figure 4. Probability of detection vs. SNR. Cooperation of 5 secondary users 
provides performance gain of 3 dB. Using multiple cyclic frequencies further 
improves the detection performance. The sum detector of the test statistic 
T> a k has the best performance. 
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Figure 5. Probability of detection vs. false alarm rate. Cooperation among 
secondary users combined with the use of multicycle sum test statistic T> a ,K 
provides the best performance. 



secondary users. Each secondary user receives the same data 
with different noise. SNR is the same for each secondary user. 

Fig. |4] depicts the performance for 5 secondary users 
compared to the single secondary user case. Performance 
gain of roughly 3 dB is obtained from the cooperation of 5 
secondary users. Using two cyclic frequencies provides similar 
performance improvement as in single secondary user case. 
Fig. [5] shows the probability of detection vs. false alarm rate 
for SNR of -9 dB. 

In the following simplistic example, we illustrate the gains 
that may be achieved via collaborative detection in the face of 
shadowing effects. In order to simulate shadowing, the SNR of 
each user was independently selected randomly from a normal 
distribution with a mean of -9 dB and standard deviation of 
10 dB. That is, the logarithm of the received power level is 
normally distributed. Fig. [6] depicts the performance of the 
multicycle detectors for the simple shadowing scenario. Com- 
parison to Fig. [5] reveals that cooperation among secondary 
users reduces sensitivity to shadowing effects significantly. 
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Figure 6. Probability of detection vs. false alarm rate. In order to simulate 
shadowing the SNR of each user was independently selected randomly from 
a normal distribution with a mean of -9 dB and standard deviation of 10 dB. 
Cooperation among secondary users reduces sensitivity to shadowing effects. 



VI. Conclusion 

In this paper, a generalized likelihood ratio test for detecting 
primary transmissions with multiple cyclic frequencies has 
been proposed, and the asymptotic distribution of the test 
statistic has been derived. In this test, impairments such as 
shadowing and fading are mitigated by combining the quan- 
tized local likelihood ratios from a number of secondary users 
under a conditional independence assumption. Simulation ex- 
amples demonstrating the improved reliability in the detector 
performance in the low SNR regime as well as significant 
gains obtained via collaborative decision making have also 
been presented. 
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